Nonlinear magnetotransport in interacting chiral nanotubes 
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Nonlinear transport through interacting single- wall nanotubes containing a few impurities is stud- 
ied theoretically. Extending the Luttinger liquid theory to incorporate trigonal warping and chirality 
effects, we derive the current contribution 1^ even in the applied voltage V and odd in an orbital 
magnetic field B, which is non-zero only for chiral tubes and in the presence of interactions. 
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Electronic transport experiments on individual single- 
wall carbon nanotubes (SWNTs) Q have revealed am- 
ple evidence for the Luttinger liquid (LL) phase of one- 
dimensional (ID) interacting metals induced by electron- 
electron (e-e) interactions. In its simplest form, the effec- 
tive low-energy theory for interacting SWNTs is in- 
sensitive to the chiral angle 0, Q describing the wrap- 
ping of the graphene sheet. This fact can be rational- 
ized by noting that, to lowest order in a k • p scheme, 
the graphene dispersion reflects an isotropic Dirac cone 
around each K point in the first Brillouin zone Q- Im- 
posing periodic boundary conditions around the SWNT 
circumference slices this cone and gives identical disper- 
sion for all 0, as long as the SWNT stays metallic. Such 
an approach is however insufficient for a description of the 
magnetotransport effects in chiral tubes. Therefore, we 
extend the theory |2| to include chirality effects by tak- 
ing into account trigonal warping, tube curvature, and 
magnetic field B, and then compute the nonlinear two- 
terminal magneto conductance. While the well-known 
Onsager symmetry G{B) = G{—B) Q excludes linear-in- 
B terms in the linear conductance, such terms can appear 
out of equilibrium [^0,I3j with first experimental obser- 
vations reported for SWNTs [lol and semiconductor 
quantum dots or rings [Til 111 lialll ITst . 

The current contribution odd in B and even in 
the voltage V is of fundamental and unique importance, 
mainly due to two reasons. First, it requires a non- 
centrosymmetric (chiral) medium, with the sign of le de- 
pending on the handedness (enantioselcctivity) [H, since 
the current density is a polar vector but magnetic field 
an axial one. Thus Je ^ requires the simultaneous 
breaking of time reversal symmetry (by the magnetic 
field) and of inversion symmetry (by the chiral medium) . 
Second; standard arguments based on the Landauer- 
Biittiker scattering formalism valid in the noninteracting 
case 0, show that 1^ ^ also requires interactions. 
At low temperature (T), e-e interactions should there- 
fore contribute to 1^ in leading order. Measurements 
of le probe interactions and chirality in a very direct 
manner, potentially allowing for the structural charac- 
terization of chiral (sin 69 ^ 0) interacting nanotubes via 
transport experiments. Nonetheless, apart from the clas- 
sical phonon-dominated high-T regime , no predictions 



specific to SWNTs have been made so far. Here we de- 
termine the current contribution (V, T) linear in the 
parallel orbital magnetic field B (Zeeman fields play no 
role here), for SWNTs in good contact to external leads. 
The simplest case allowing for Jg 7^ is found when in- 
cluding at least two (weak) elastic scatterers, represent- 
ing either defects or residual backscattering induced by 
the two contacts. Our theory includes the often strong 
e-e interactions using the bosonization method 0|, and 
holds for arbitrary chiral angle 9. Interactions are char- 
acterized by the LL parameter if < 1 in the charge sec- 
tor, see Eq. (j^Jl below, with typical estimate K « 0.2 
Our result for le, see Eqs. @) and H12|) . is non-zero 
only for chiral (sin 661 7^ 0) interacting {K < 1) SWNTs, 
and changes sign for different handedness {6 —6'). We 
predict oscillatory behavior of Ig as a function of bias 
voltage, where the oscillation period depends on the Lut- 
tinger parameter K, see Eq. (|15|l . The coefficient ^3 
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is shown to exhibit power-law scaling at T ^ 0, with a 
negative exponent that is again determined by K. In- 
cluding only elastic impurity backscattering, a(T) does 
not change sign with T for otherwise fixed parameters. 

Chirality effects in the low-energy theory come about 
when one includes trigonal warping in the band struc- 
ture. Its main effect is to introduce different Fermi ve- 
locities Vra |l7l | for right- and left-moving excitations 
(r = R/L = ±) around the two distinct K points 
(a ~ ±). All these velocities coincide when disregard- 
ing trigonal warping, and in a nominally metallic SWNT 
are given by Vra — v = 8x 10^ m/sec. For clarity, we con- 
sider an electron-doped SWNT with equilibrium chemi- 
cal potential fi = hvkp > well inside the conduction 
band, and omit the valence band. Including the trigonal 
warping to lowest nontrivial order U 0] , we find 
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The effects of tube curvature and magnetic field enter via 
the quantized transverse momenta 



k^aR = n- av/i + $/$o + accos(36'), 



(3) 



where a = 0.246 nm is the lattice spacing, R the tube 
radius, <I> = ttR^B, <i>o = h/e the flux quantum, and the 
tube curvature results in c = na/R with k « 1 Here 
integer n with \k^a\ < kp are allowed; we focus on the 
lowest band n = in what follows. The index v distin- 
guishes nominally metallic {v = 0) and semiconducting 
{v = ±1) SWNTs, and for simplicity, from now on we 
assume v = Q. However, with minor modifications, our 
theory below also applies to strongly doped semiconduct- 
ing tubes. We mention in passing that chirality effects 
are also important for other physical quantities. In par- 
ticular, scattering by a long- range p otential in metallic 
SWNTs depends on chirality ^t^. 

In terms of annihilation fermion operators Raa{x) and 
Laa{x) for right- and left-movers of spin a = ±, re- 
spectively, the usual linearization of the band structure 
around the Fermi points then yields H = iJ^L + Vdis (we 
set h = ks = 1), 

HlL = {""Ra^iadxRaa - '"LaLi^dxL^^) 

+ Y j dx{R^R + L^L)\ (4) 
Vd^s = / C/(a;) (e-2''=^"i?tL + h.c.) , (5) 



where 4Vb /ttv = — 1 describes e-e forward scattering 
interactions pol |. and aa summations are implied when 
not given explicitly. Elastic disorder leads to V^is, e.g., 
due to impurities, defects or substrate inhomogeneities. 
We keep only the intra-band backscattering potential 
C/(x), which yields the dominant impurity effect 0. 
This Hamiltonian can be efficiently treated by (Abelian) 
bosonization |l6j |. Introducing four bosonic fields (j)i{x) 
and their dual 0i{x), with i = (c+, c—, s+, s—) denoting 
the total/relative charge/spin modes, the clean Hamilto- 



nian for Vr, 
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J2 J dx + (8,0,)'] , (6) 
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where Uc = v/ K is the plasmon velocity for the c-|- mode, 
and Un — V for the three neutral modes. Including the 
Vra differences in Eq. ((5J then brings about two new fea- 
tures, acting separately in the decoupled charge {i ~ c±) 
and spin (i = s±) sectors: (i) couplings between total 
{+) and relative (— ) modes, and (ii) couplings between 
mutually dual {9i, <f>i) fields. The first leads to tiny quan- 
titative corrections but no qualitative changes and is ne- 
glected henceforth. Point (ii) is crucial, however, since it 
implies different velocities {ucr/l a-nd Un^R/i,) for right- 



and left-moving plasmons. This eventually produces the 
current contribution 1^ in the presence of impurities. To 
linear order in i?, we find 
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Uum/l/v = '^±5, (7) 
[a/Rf sin(66'). 



Note that (5 7^ requires both B ^ and chirality, 
sin 60 7^ 0. Moreover, 5 has opposite sign for oppo- 
site handedness. It essentially describes the difference 
of i?/i-moving velocities, and depends linearly on B. 

Next we address the computation of Je(y,T), where 
the current operator is {2e / y^)dt4>c+ ■ An important 
point concerns the inclusion of the applied voltage V in 
a two-terminal setup. For weak impurity backscattering 
U{x), it is sufficient to address the clean case with adia- 
batically connected leads, where a time-dependent shift 
arises |2j, (pc^ — > + eVt/y/r:. However, with chiral 
asymmetry, vr 7^ vl, there is an additional effect [23 | 
due to the different density of states i'r/l = 2/7rUfl/^ 
for R/L movers. Starting from the equilibrium chemi- 
cal potential /.(, when applying a voltage, the chemical 
potentials fJ.R/L of R/L movers are set by the left /right 
reservoirs, respectively, where fiR ~ — sV . Relative 
to the equilibrium density, R/L movers arc thereby in- 
jected with density p^fyi ~ ^r/ L{^J^R/ l ^ m)- We may 
now write ^jlr/l — + ^ 6^/2, where A/i is de- 
termined by the condition that in an adiabatically con- 
nected impurity- free quantum wire, no charge can accu- 
mulate in the steady state |23|, Pr+ Pl = 0- This yields 
A/i = SeV/2, where 5 = {vr — vl)/{vr + vl) has been 
specified in Eq. |(7J), implying the voltage-dependent shift 
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kF + 5eV/2vo, vo = ^^^^^ ^ v{l - S^) (8) 
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in Eq. ((SJ. Under the non-equilibrium Keldysh formal- 
ism, to second order in U (x) , 1^ then follows as 

2eK f°° f 
le = / dt / dXdxsm{2kFx) (9) 

[nay J_oo J 

X U{X + x/2)U{X -x/2) 

X cos(ey[t + ,5x/uo]) Ime*'^^^^^'*), 

where is the sum of the lesser Green's functions for 
the four plasmonic modes, see Eq. 1)11(1 below. 

As a function of time, e'-'^S^{x,t) j^^^g singularities only 
in the lower-half complex plane, and therefore /g = for 
V — 0. By inspection of Eq. we also observe that 
(5) = With Eq. (O we conclude that Ig is 

odd both in magnetic field and chiral angle. In fact, le is 
the only odd-in-i? contribution, since the odd-in-l^ part 
of the current turns out to be even in B. We can now 
state a first necessary condition for 1^ 0, namely (5 7^ 0, 
which requires the simultaneous breaking of time reversal 
and inversion symmetry. Approximating the impurity 
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potential as U{x) = J^i Ui5{x — Xi), Eq. © becomes 
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with Xjk = Xj — Xk- Then a second necessary condi- 
tion arises: there should be just a few impurities. The 
function Q^{x,t) changes slowly on the scale kp^, and 
when one has to sum over many impurities, 1^ = due 
to the fast oscillations of s'm{2kpXjk) in Eq. (|10|) . In 
such a disordered case, however, Ructuations exhibit re- 
lated magnetochiral transport effects Q. When only a 
few impurities (but at least two) are present, there is no 
averaging and the effect survives in the current, albeit its 
magnitude and sign are of course sample-dependent. It 
is likely that the experiments of Ref. [13 were performed 
on samples with not too many impurities, for otherwise 
strong localization effects characteristic for ID systems 
would render them insulating. Below we focus on the 
simplest case of two impurities separated by a distance 
L ~ X2 — xi, see also Ref. j^^. Such a controlled two- 
impurity setup can be experimentally realized in individ- 
ual SWNTs Ui- 

As further evaluation of Eq. H1U|) requires numerical 
integration routines, we first examine the spinlcss single- 
band version of Eq. Q), which contains the essential 
physics and allows to compute /g in closed form |25| . 



Given the Fermi velocities Vr- 



R/L=± 



V{1 ± S) 



corre- 



sponding to Eq. (|2J)] , the plasmon velocities of the clean 
interacting system are Ui?/L = v{l/K±6) [corresponding 
to Eq. ((TJ]. The plasmon lesser Green's function is 
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Using 2uiiul/ {ur + Up) ~ v/K, some algebra yields 
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where V is the Gamma function and F the hypergeomet- 
ric function p^ . We introduced dimensionless tempera- 
ture and voltage 
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Obviously, /g = in the noninteracting case {K — 1), 
and thus interactions provide a third necessary condi- 
tion for /e 7^ 0. Equation l|12|) now allows to analyze 




FIG. 1: Even-in-voltage part of the current le (in nA) as 
a function of bias voltage (in V) for several temperatures. 
Results are shown for a (10,4) SWNT with K = 0.23, B = 
16 T, containing two impurities (L = 20 nm, f/1,2 = t'/2). 
Inset; Coefficient a{T) [in ^SV^T''^] defined in Eq. as 
a function of temperature in double-logarithmic scales. 



several limits of interest. First, one recognizes an oscilla- 
tory dependence on the doping level ^ = vkp tuned by a 
backgate voltage, similar to what is seen experimentally 
[Toj . Second, a{T) [see Eq. (0)] is exponentially small for 
T 1 but shows power-law scaling a{T <C 1) oc T'^^-^ 
at low temperatures, implying a huge increase in a{T) 
when lowering T. Third, in the zero-temperature limit, 
Eq. lO yields 

Ie{V) cx sin[(5(l - K^)V/K] V""-'/^ Jk-i/2{V) (14) 

with the Bessel function Ji, j^^. This implies the low- 
voltage scaling If, cx V^^ , a posteriori justifying our per- 
turbative treatment of the impurity potential. Remark- 
ably, /e shows oscillatory behavior as a function of the 
bias voltage V. As one can see from Eq. (|14|l . there are 
two different oscillation periods. 
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For strong interactions, AV2 can in principle provide di- 
rect information about S. In any case, AVi should be 
readily observable and yields already the Luttinger pa- 
rameter K. Inspection of Eq. ((Tn|l suggests that the phys- 
ical reason for these oscillations is the quantum inter- 
ference between right- and left-moving waves travelling 
between the two impurities with different velocities. 

Let us then go back to the full four-channel case, and 
perform the integration in Eq. H1Q(I numerically. To be 
specific, we consider a long and adiabatically connected 
(10,4) SWNT with two symmetric impurities separated 
by L = 20 nm. Fig.^shows /g as a function of V for sev- 
eral T. Oscillatory behaviors with bias voltage are clearly 
visible. With decreasing temperature, Ig increases and 
one gets power-law scaling of a{T) (x T^^"""^^/^ at low T 
(see inset of Fig. generalizing the above single-channel 
result. Notably, a{T) has the same sign for a given pa- 
rameter set, in qualitative agreement with experiments 
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[lOj- Within our parameter choiees, also the order of 
magnitude in a{T) agrees with Ref. case, 
the picture obtained in the single-channel version is es- 
sentially recovered under the four-channel calculation. 

To conclude, we have analyzed nonlinear magnctochi- 
ral transport properties of interacting single-wall carbon 
nanotubcs. In chiral tubes, measurement of the odd-in-_B 
component Ig (which must be even in V) provides direct 
information about interactions and chirality not acces- 
sible otherwise. For two impurities, we have presented 
detailed analytical results for If.- We predict oscillations 
of le as a function of bias voltage, which provide direct 
information about the interaction parameter K . More- 



over, at low temperatures, power-law scaling is found and 
leads to an enhancement of Ig- In future work, our ap- 
proach should be useful when calculating fluctuations of 
transport coefficients in disordered interacting SWNTs. 
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